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^ ■ Abstract— Currently available data on the field of velocities VJ, H for open star 

clusters are used to perform a kinematic analysis of various samples that differ by heliocentric 
^ . distance, age, and membership in individual structures (the Orion, Carina-Sagittarius, and 

Perseus arms). Based on 375 clusters located within 5 kpc of the Sun with ages up to 



> 



O I 1 Gyr, we have determined the Galactic rotation parameters uq = —26.0 ±0.3 km s~^ kpc^^, 

■ i^o ~ 4.18 ± 0.17 km s~^ kpc~^, uJq = —0.45 ± 0.06 km s~^ kpc~^, the system contraction 

parameter K = —2.4 ± 0.1 km s^^ kpc^^, and the parameters of the kinematic center 
Rq = 7.4 ± 0.3 kpc and Iq = ± 1°. The Galactocentric distance Rq in the model used has 
been found to depend significantly on the sample age. Thus, for example, it is 9.5 ± 0.7 
kpc and 5.6 ± 0.3 kpc for the samples of young (< 50 Myr) and old (> 50 Myr) clusters, 
respectively. Our study of the kinematics of young open star clusters in various spiral arms 
has shown that the kinematic parameters are similar to the parameters obtained from the 
Q\ • entire sample for the Carina-Sagittarius and Perseus arms and differ significantly from them 

O I for the Orion arm. The contraction effect is shown to be typical of star clusters with various 

Q I ages. It is most pronounced for clusters with a mean age of ~ 100 Myr, with the contraction 

velocity being Kr = —4.3 ± 1.0 km s~^. 

S INTRODUCTION 

The Galactic rotation parameters have been repeatedly determined by many authors using 
objects belonging to various Galactic structural components: from ionized and neutral hy- 
drogen (Fich et al. 1989; Merrifield 1992; Brand and Blitz 1993; Nikiforov 1999; Avedisova 
2005), from distant OB associations of stars (Dambis et al. 2001; Mel'nik et al. 2001), and 
from open star clusters (Zabolotskikh et al. 2002; Gerasimenko 2004; Popova and Loktin 
2005b). 

Open stars clusters (OSCs) are of great interest in studying the kinematics of the Galaxy, 
since they are located in a wide solar neighborhood and have reliable distance and age 
estimates. 

Data only on the radial velocities Vr (Mishurov et al. 1997; Gerasimenko 2004; Popova 
and Loktin 2005b) are commonly used for a kinematic analysis of OSCs, because the random 
errors in the radial velocities are essentially distance-independent and, hence, distant objects 
can be used. Note also that data only on the radial velocities are available for distant 
hydrogen clouds (Fich et al. 1989; Avedisova 2005). 
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When such catalogs as Hipparcos (1997) and Tycho-2 (Hog et al. 2000) appeared, it 
became possible to accurately determine the mean proper motions of OSCs (Beshenov and 
Loktin 2004; Kharchenko et al. 2005b). A number of authors used simultaneously two 
observed velocity components, Vr and Vi (Mishurov and Zenina 1999; Dambis et al. 2001), 
or all three velocity components, Vr, Vi, and Vb (Zabolotskikh et al. 2002; Bobylev 2004), 
to samples of stars within 3.4 kpc of the Sun. In this paper, we also used all three observed 
velocity components. 

One of the important problems is to determine the distance to the center of Galactic 
rotation Rq. This parameter is estimated indirectly from an analysis of the velocities of 
objects by reconciling the adopted model with observational data. Only in recent years has 
this parameter been estimated directly from objects of the Galactic nucleus (McNamara 
2000; Eisenhauer et al. 2003). A summary of the Rq values determined in the last decade 
by various methods can be found in Avedisova (2005) and an overview of the Rq estimation 
methods is given in Nikiforov (2004). 

The goal of this paper is to determine the Galactic rotation parameters {uq, uj'q, Uq) and 
the parameters of the kinematic center (Galactocentric distance Rq and direction Iq) from 
currently available data on the field of space velocities Vr, V, and Vb, distances, and ages of 
OSCs and to study the kinematic peculiarities of various samples differing by heliocentric 
distance, age, and membership in individual arms. 

DATA 

At present, more than 1700 OSCs are known in the solar neighborhood. Data on their 
proper motions, radial velocities, and positions are needed for our purposes. A catalog that 
includes 652 OSCs (Kharchenko 2001; Kharchenko et al. 2005a, 2005b; Piskunov et al. 2006) 
forms the basis for our work list. The advantage of this catalog is a homogeneity and a high 
accuracy of the determination of mean cluster proper motions achieved by using the ASCC- 
2.5 all sky catalog (Kharchenko 2001) compiled from Hipparcos (1997), Tycho-2 (Hog et al. 
2000), and several other sources. The age estimates obtained by comparison with isochrones 
are available for the clusters of this catalog. The cluster distance estimates are based on the 
results by Loktin and Beshenov (2003), who reconciled the photometric estimates with the 
Hipparcos distance scale. We took other data from the compilation by Dias et al. (2002) and 
the WEBDA database (http://obsww.unige.ch/webda/). 

For such open clusters as ASCC 16, ASCC 18, and Tr 10, we used the mean radial 
velocities that we improved (Bobylev 2006) using the OSACA catalog (Bobylev et al. 2006a). 

As a result, we compiled a database on the proper motions, radial velocities, and distances 
of 394 OSCs. They are located within about 5 kpc of the Sun. Their ages do not exceed 
1.5 Gyr. The radial velocities with estimates of their random errors are available for 270 
(~ 70% of the sample) clusters. Note for comparison, that Dias and Lepine (2005) used only 
212 OSCs with measured proper motions and radial velocities served as a source for their 
kinematic analysis. 

To identify clusters belonging to various structures (e.g., the Orion, Perseus, and Carina- 
Sagittarius arms), we use a probabilistic approach that we developed (Bobylev and Bajkova 
2007). It is based on the approximation of the two-dimensional {XY) probability density 
function for the objects under consideration by a set of Gaussians that represent the proba- 
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bility density functions of individual features. 



THE METHOD 

In this paper, we use a rectangular Galactic coordinate system with the axes directed away 

from the observer toward the Galactic center {I = 0°, b = 0°, the X axis), along the Galactic 
rotation (/ = 90°, b = 0°, the Y axis), and toward the North Galactic Pole {b = 90°, the Z 
axis) . 

The youngest disk objects that lie in the Galactic plane are commonly used to determine 
the Galactic rotation parameters. As a result, the observed velocity component H is generally 
not considered, while the component of solar motion Wq is assumed to be known (Dambis 
et al. 2001; Gerasimenko 2004). At the same time, analysis of relatively distant objects 
(Zabolotskikh et al. 2002; Bobylev 2004) showed that including the components Wq and H 
in the model allows their values to be also determined rehably. Therefore, we dwell on this 
more general approach. In the special case where the clusters are distributed in a ring (see 
the Carina-Sagittarius and Perseus sample below), we used a fixed value of Wq — 7.2 km s"^ 
(Dehnen and Binney 1998). 

The method for determining the kinematic parameters used here consists in minimizing 
the quadratic functional F, 

N N N 

min F = - v:f + E^iW - v;r + E^^(K - vif (1) 

1=1 i=l i=l 

under the following constraints derived from Bottlinger's formulas (Ogorodnikov 1965) with 
the angular velocity of Galactic rotation expanded in a series to terms of the second order 
of smallness in t/Rq: 

Vr — —Uq cos6cos(Z — Iq) — (2) 
—Vq cos b sin(Z — Iq) — Wq sin b— 
—Ro{R — Rq) sm{l — Iq) cosbujQ— 
-0.5i?o(i? - Ro)^ sin(Z - Iq) cos 64'+ 
+rK cos^ b, 

Vi^Uq sin(Z - Iq) - Vq cos(Z - Iq) - (3) 
— {R — Ro){Ro cos(/ — Iq) — r cosb)ujQ— 
-{R- Ro)^{Ro cos{l - Iq) -rcosb)x 
x0.5cJq + rujQ cos b, 
Vb = Uq cos(/ - /o) sin 6 + (4) 
+Vq sin(/ — Iq) sin b — Wq cos b+ 
+Rq{R — Rq) sin(/ — Iq) sin 6cJq+ 
+0.5Rq{R - Rq)^ sin(/ - Iq) sin bujQ- 
—rK cosbsinb, 

where is the number of clusters used; i is the current cluster number; are the cluster 
velocities to be calculated, with Vr being the radial velocity and Vi — 4.74r//;COs6 and Vj, = 



A.lAjjLb being the proper motion velocity components in the / and h directions, respectively 
(the coefRcient 4.74 is the quotient of the number of kilometers in an astronomical unit by the 
number of seconds in a tropical year); V^^Vl, are the measured components of the velocity 
field; Wr,wi,Wh are the weights; r is the heliocentric distance of the cluster calculated via the 
photometric parallax determined by bringing the cluster main sequence into coincidence with 
the corresponding isochrone (Piskunov et al. 2006); the cluster proper motion components 
III cosfo and fj,b are in mas yr~^, the radial velocity Vr is in km s~^; Uq, Vq, Wq are the cluster 
centroid velocity components relative to the Sun; Rq is the distance from the Sun to the 
kinematic center of the system; R is the distance from the cluster to the center of rotation; 
Iq is the direction of the kinematic center; R, Rq, and r are in kpc. The quantity ujq is 
the angular velocity of rotation at distance Rq, the parameters uj'q and lUq are the first and 
second-order derivatives of the angular velocity, respectively; K is the Oort constant that 
describes the expansion/contraction of the stellar system. The distance R can be calculated 
using the expression 

i?2 ^ cos bf - 2Ror cos b cos(/ - Iq) + Rl- (5) 

Note also that Eqs. (2)-(4) are written in such a way that the direction of rotation from the 
X axis to the Y axis is positive. 

The weights in functional (1) are assigned in accordance with the expressions (for sim- 
plification, the index i was omitted) 

= VA' + ^v., (6) 

where 5*0 denotes the dispersion averaged over all observations, which has the meaning of 
the "cosmic" dispersion that wc take to be 8 km s~^; f] = av^/cry^ and 7 = ay^/av^ are the 
scale factors (in our case, j3 = 1 and 7 = 2). The errors in the velocities Vi and VJ, can be 
calculated using the formula 

where aT^/ir is taken to be 0.2, a typical error in the photometric parallax. 

In this paper, apart from the system of weights (6) described above, we also use a variant 
of unit weights where Wj. = Wi = wj, = 1 for comparison. 

The optimization problem (l)-(5) is solved for nine unknown parameters Uq, Vq, Wq, luq, 
lu'q, uJq, K, Rq, and Iq by a coordinate- wise descent method (the sought-for parameters are 
taken as the coordinates). According to this method, the increment Zx for each unknown, 
which we arbitrarily denote by x, is sought from a necessary condition for the existence of 
an extremum of the functional Flzx). 

^ = 0. (7) 

dzx 

The functional F{zx) is obtained by substituting Eqs. (2)-(4) for Vr{zx), Vi{zx), Vb{zx) into 
(1) using (5) in which the variable x is replaced hy x + Zx- 
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Generally, Eq. (7) is a nonlinear equation for Zx- Its numerical solution by the Newton 
method is an iterative process. The value for Zx at iteration i is 

4 = - {dF{zx)/dZx)/{d?F{zx)/dz'^ I 

The solution for Zx usually converges after 1-3 iterations. As a result of one descent, the 
sought- for variable acquires a new value, x + Zx- In our case, about 1000 descents in all nine 
sought-for variables were required to obtain a solution with a sufficiently high accuracy. 

A sufficient condition for the existence of a global extrcmum is that the Hessian ma- 
trix composed of elements {aij} = d'^F/dxidxj, where Xi{i = 1, ...9) denote the sought-for 
parameters, be positively defined everywhere. We calculated the Hessian matrix in a wide 
domain of parameters: 4 < Uq, Vq < 16, 2 < < 11, —30 < a;o < 0, — 3 < a;Q < 6, —0.6 < 
cUq < 0.3, —4 < K < 5,3 < Ro < 15, —20 < lo < 10. Our analysis of the Hessian matrix 
for both weighting variants for all of the cluster samples considered here showed it to be 
positively defined. This suggests the existence of a global minimum in this domain and, as a 
result, the uniqueness of the solution. For unit weights, the Hessian matrix is also positively 
defined far outside this domain, which could not be established in the case of weighting 
according to rule (6). However, as will be shown below, the adopted weighting allowed the 
accuracy of the solutions obtained to be increased noticeably. 

Interestingly, using weights of the form Wr,i,b = l/o'y provided no positive definiteness 
of the Hessian matrix in the domain considered. It thus follows that several local extrema 
exist in this case. Therefore, to find the correct solution, we must know the initial approx- 
imation as accurately as possible. In contrast, if a global extremum exists, then we can 
start from any point of the domain. This, in particular, explains the modification of the 
weights according to rule (6), which provided a global minimum in a fairly wide domain of 
parameters. 

We estimated the errors in the sought-for parameters by means of Monte-Carlo simula- 
tions. The errors were estimated by performing 100 cycles of computations. At this number 
of cycles, the mean values of the solutions virtually coincide with the solutions obtained only 
from the input data, i.e., without any addition of the measurement errors. 

CONSTRAINTS 

In solving the optimization problem (l)-(5), we used the following constraints on the data: 
(i) the magnitude of the pecuhar velocity \Vpec\ < 90 km s~^; (ii) the cluster height above 
the Galactic plane \z\ < 500 pc. The velocity Vpec is calculated with respect to the local 
standard of rest -^0 = 10.0 km s~^, Vq = 5.3 km s^^, Wq = 7.2 km s^^ (Dehnen and Binney 
1998). According to constraint (i), four clusters were rejected. Constraint (ii) is important 
only for clusters older than 50 Myr. According to this criterion, the seven oldest clusters 
were rejected; all of the remaining clusters proved to be younger than 1 Gyr. After applying 
an additional constraint on the sample radius, r < 5 kpc, 375 clusters remained out the 
originally compiled list of 392 clusters, which served as the material for a further analysis. 
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RESULTS 



Figure 1 presents the spatial distribution of OSCs in the Galactic XY plane. Since a kine- 
matic analysis of OSC samples with age separation is also the goal of this paper, we separately 
displayed the distributions of OSCs younger and older than 50 Myr in Figs, la and lb, re- 
spectively. We see from Fig. la that young clusters clearly reproduce three segments of the 
spiral arms (Mel'nik et al. 2001; Dias and Lepine 2005; Popova and Loktin 2005a). The 
concentration regions are marked in the figure by the ellipses: the Perseus arm is on the left, 
the Orion arm is at the center, and the Carina-Sagittarius arm is on the right. We clearly see 
that the distribution of older clusters (Fig. lb) is considerably more uniform, more compact, 
and concentrated closer to the Sun. 

Figure 2 shows the mean errors cry^ , , and ay^ as a function of the heliocentric distance. 
Note that in solving the optimization problem (l)-(5), we take into account the peculiarities 
of the error distribution using the weighting procedure (6). 

Table 1 presents the results of our kinematic analysis of OSC samples with different 
sample radii (r < 5 kpc and r < 2.5 kpc) and age separation (younger and older than 50 
Myr) using the two weighting methods. Column 1 gives the sample type and size; column 
2 lists the centroid velocity components relative to the Sun; columns 3, 4, and 5 list the 
rotation parameters; columns 6, 7, and 8 give the values of the ii'-effect, Ro, and Iq. The 
results obtained with unit weights and weights (6) are presented in the upper and lower parts 
of the table, respectively. As can be seen from Table 1, the random errors in virtually all 
parameters decrease when the system of weights (6) is used. For both weighting methods, 
the difference in Rq, which is significant for the sample of old clusters, is largest. As can 
be seen from Table 1, we obtain a discrepancy in Uq depending on the constraint on the 
sample radius. Figure 3 shows how the parameter cUq affects the Galactic rotation curve. To 
properly compare the results and to construct the curves indicated in Fig. 3, we obtained 
two solutions at fixed Rq — YA kpc: 

cuo = -26.7 ± 0.3 km s'^ kpc~\ (8) 
lu'q = 4.27 ±0.11 km s^^ kpc'^ 
Uq = -1.05 ± 0.06 km s'^ kpc'^, 

for the r < 2.5 kpc sample and 

= -26.0 ± 0.2 km s'^ kpc-^ 
uj'q = 4.14 ± 0.10 km s^^ kpc'^ 
oj'^ = -0.44 ± 0.04 km s"^ kpc'^ 

for the r < 5 kpc sample. As can be seen from Fig. 3, the difference in the rotation curve 
due to the influence of the second derivative becomes noticeable as one receded from Rq to 
distances exceeding pa 1.5 kpc. 

The Oort constants A — O.SRqUJq and B — u;o + O.^RqUq calculated from the data in the 
fifth row of Table 1 are >1 = 15.4 ± 0.6 km s"^ kpc"^ and B = -10.6 ± 0.7 km s"^ kpc"^ 
and have no significant differences for all solutions of Table 1. 

In Fig. 4, linear contraction velocity Kr is plotted against heliocentric distance r and 
mean sample age t. To obtain these characteristics, we divided the entire data set into four 
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groups in such a way that each group contained approximately the same number of clusters. 
The dependences in Fig. 4 were derived at fixed Galactocentric distance Rq — 7.4 kpc. It 
follows from Fig. 4a that the contraction effect is traceable at all hehocentric distances r > 1 
kpc and is, on average, Kr = —3.2 ± 1.0 km s~^. As we see from Fig. 4b, the contraction 
effect is typical of star clusters with various ages, but it is most pronounced for clusters with 
a mean age t ~ 100 Myr; the contraction velocity in this case is Kr = —4.3 ± 1.0 km s~^. 

To study the kinematic peculiarities of various structures, we analyzed clusters younger 
than 50 Myr belonging to various arms (Fig. la). To identify the arms, we calculated the 
probability that each cluster belonged to one or another arm and then attributed it to the 
arm the probability of belonging to which was at a maximum. In this case, the probability 
densities of cluster belonging to each arm were fitted by Gaussians whose bases in the shape 
of ellipses are shown in Fig. lb. The results of our kinematic analysis in which all parameters 
were found simultaneously are presented in the first four rows of Table 2. The first and second 
two rows give, respectively, the results obtained with unit weights and with weights (6). 

For the Carina-Sagittarius and Perseus sample, which contains mostly distant clusters, 
we also obtained a solution at fixed Wq = 7.2 km s~^ using weighting (6). This solution is 
given in the last row of Table 2. It shows that a change in parameter Wq by ~ 2 km s~^ 
does not affect significantly the determination of other model parameters. 

DISCUSSION 

Parameters of the Kinematic Center 

As can be seen from Table 1, the values of Rq obtained with various weighting methods 
agree well, within the la error limits. Rq — 7.4 ± 0.3 kpc that wc found using clusters with 
a mixed (in age) composition is in good agreement with the Rq estimates made by various 
authors. 

Thus, for example, in the opinion of Avedisova (2005), the most probable value of Rq 
lies within the range from 7.5 to 8.2 kpc. By analyzing the Rq determinations by various 
authors, Nikiforov (2004) showed that the currently best vahic of Rq is 7.9 ± 0.2 kpc. Using 
270 OSCs, Chen and Zhu (2007) obtained an estimate of Rq = 7.95 ± 0.62 kpc. 

Note that the Galactocentric distance Rq = 7.4±0.3 kpc we found provides more evidence 
for the short distance scale (Dambis et al. 2001) than the lAU recommendations (1985), 
Rq — 8.5 kpc. 

At the same time, as we see from Table. 1, Rq depends significantly on the sample age. 
Thus, for example, it is Rq = 9.5 ± 0.7 kpc for clusters younger than 50 Myr. On the other 
hand, a value of -Rq = 7.1 ± 0.6 kpc, which is close to that found from all 375 OSCs, was 
found from clusters younger than 50 Myr that belong only to the Carina-Sagittarius and 
Perseus arms, but not to the Orion arm (Table 2). Our separate study of the Orion arm 
showed a significant difference almost in all kinematic parameters (see Table 2). 

Clusters older than 50 Myr are distributed fairly uniformly in space (Fig. lb); the 
influence of the Galactic spiral structure and the Local system of stars in their motions is 
minor. One would think that they should be well suited to a reliable determination of Rq, 
but we found Rq — 5.6 ± 0.3 kpc for them (the lower part of Table 1). Based on a kinematic 
analysis of the OSC radial velocities from various catalogs and, in particular, the catalog by 
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Piskunov et al. (2006), Nikiforov and Kazakevich (2006) point out a great uncertainty in 
Rq. In particular, they found i?o = 6 ± 0.7 kpc from a sample of old (logi > 8.8) clusters, 
in good agreement with our result. 

Analysis of the motions of the stars nearest the Sun showed that the vertex deviation 
for the youngest stars reaches 30° (Dehnen and Binney 1998). Based on the clusters the 
results of whose analysis are presented in Table 1, we found no significant deviation from the 
direction Iq — 0°. This means that the rotation of the Galactic disk is nearly axisymmetric 
on large scales. 

The Galactic Rotation Curve 

The parameters of the Galactic rotation curve that we found in solution (8) at fixed Rq — 7.4 
kpc (which, in our case, is best suited to a proper comparison) are in good agreement with 
coo = -27.5 ± 1.4 km s'^ kpc-\ cu'^ = 4.54 ± 0.24 km s'^ kpc'^, u'^ = -1.09 ± 0.19 km s'^ 
kpc~^ determined from open clusters by Zabolotskikh et al. (2002) for the short distance 
scale and for Rq = 7.5 kpc. There is also good agreement with the results of our analysis of 
distant OB stars (Bobylev 2004): = -28.0 ± 0.6 km s^^ kpc-\ uj'o = 4.17 ± 0.14 km s"^ 
kpc-2, u'^ = -0.81 ± 0.12 km s-^ kpc'^ (for Ro = 7.1 kpc). 

We found the parameter cJq to decrease in absolute value with increasing sample radius. 
Far from the Sun, the Galactic rotation curve constructed from our data (Fig. 3) disagrees 
with the rotation curve constructed, for example, from the molecular gas (Avedisova 2005), 
which is flat in a wide R range, from 2 kpc to 15 kpc. Therefore, using a larger number of 
terms in the expansion of cuq in terms of r/Ro is of considerable interest in establishing the 
shape of the Galactic rotation curve from OSCs far from the Sun. 

The Contraction Velocity 

It is interesting to compare the contraction parameter K = —2.4 ±0.1 km s~^ kpc~^ that we 
found from all clusters with the results by Torra et al. (2000). Based on the Ogorodnikov- 
Milne linear model, these authors found K = —2.0 ± 0.4 km s^^ kpc^^ for OB stars of all 
ages in the range of distances 100-2000 pc and K = —5.1 ± 1.5 km s^^ kpc^^ for OB stars 
with ages t > 60 Myr in the range of distances 100-2000 pc. Fernandez et al. (2001) and 
Bobylev et al. (2006b) showed that including the spiral structure did not eliminate the K- 
effect. Based on our results (Fig. 4), we can assume that the observed contraction velocity 
reflects the fllling rate of the interarm space (Marochnik and Suchkov 1984), but a detailed 
study of this effect is outside the scope of our paper. 

Kinematic Peculiarities of the Arms 

In this paper, wc applied no corrections for the influence of the spiral structure to the 
observed velocities. However, the influence of the spiral structure was studied in our previous 
paper (Bobylev et al. 2006b), where we showed that it has a major effect on the centroid 
velocity components relative to the Sun. This can explain the significant difference in the 
components A-U0 ^ 2 km s~^ and Avq ~ 5 km s~^ found from various arms. 

Analysis of the kinematic parameters for young clusters belonging to the Orion arm 
indicates that such parameters as ujq = —31.7 ± 0.6 km s~^ kpc^^ and Iq = —3 ± 1° differ 
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significantly from tlie Galactic rotation parameters found from all clusters. This shows that 
the Orion arm has an additional rotation with an angular velocity of ~ 5 km kpc~^ 
around a center that may not be associated with the Galactic center. 

CONCLUSIONS 

We compiled a database on 375 OSCs from various present-day catalogs and published 
sources. The clusters considered are located with < 5 kpc of the Sun. Their ages do not 
exceed 1 Gyr. 

The following Galactic rotation parameters were determined from data on the field of 
velocities K, Vu H for these OSCs: ujq = -26.0 ±0.3 km s"^ kpc-^ co'^ = 4.18 ±0.17 km s'^ 
kpc~^, Uq = —0.45 ± 0.06 km s~^ kpc^^, the system contraction parameter K = —2.4 ±0.1 
km s~^ kpc~^, and the parameters of the center of rotation Rq = 7.4 ± 0.3 kpc and Iq = 0°. 

In addition, we performed a kinematic analysis of various OSC samples differing by 
heliocentric distance, age, and membership in individual arms. The value of Rq was found 
to depend on the sample age. Thus, for example, it is 9.5 ±0.7 kpc for clusters younger 
than 50 Myr and 5.6 ± 0.3 kpc for clusters older than 50 Myr. Our study of the kinematics 
of young OSCs in various arms showed that the derived kinematic parameters are similar to 
the parameters obtained from the entire OSC sample for the Carina-Sagittarius arms and 
differ significantly from them for the Orion arm. The contraction effect was shown to be 
typical of star clusters with various ages and to be most pronounced for clusters with a mean 
age of ~ 100 Myr; the linear contraction velocity in this case is Kr — —4.3 ± 1.0 km s~^. 
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Table 1: Kinematic parameters of Galactic rotation 



Sample 


Uq,Vq,Wq, 

km/s 


t^O, 

km/s/kpc 


^0' 

km/s/kpc^ 


km/s/kpc^ 


K, 
km/s/kpc 


Ro, 
kpc 


^0, 

deg. 


r < 5 
kpc 
(N=375) 


9.5(0.4) 
10.9(0.3) 

9.1(0.2) 


-26.7(0.4) 


3.91(0.21) 


-0.35(0.08) 


-2-3(0.1) 


8-0(0.4) 


0(1) 


r < 2.5 

kpc 

(X=34(J) 


9.2(0.3) 

11.9(0.3) 

8.1(0.2) 


-27.2(0. .5) 


3.72(0.32) 


-0.70(0.1,3) 


-2.6(0.2) 


8.5(0.7) 


^2(1) 


t < 50 
Myr 
(N=196) 


9.6(0.4) 

11.7(0.4) 

9.4(0.3) 


-27.4(0.5) 


3.21(0.25) 


-0.14(0.07) 


-1.4(0.2) 


9-8(0.7) 


0(1) 


t > 50 

Myr 
(N=178) 


9.7(0.4) 
10.2(0.4) 
8.8(0.4) 


-25.3(0.9) 


4.59(0.60) 


-0.59(0.19) 


-3-5(0.3) 


6-7(0.8) 


0(1) 


r < 5 
kpc 
(N=375) 


9.2(0.2) 
10.7(0.2) 
8.2(0.2) 


-26.0(0.3) 


4.18(0.17) 


-0.45(0.06) 


-2.4(0.1) 


7-4(0.3) 


0(1) 


r < 2.5 

kpc 
(N=340) 


8.8(0.2) 
12.0(0.2) 
8.0(0.1) 


-26.7(0.3) 


3.45(0.26) 


-0.69(0.11) 


-2-4(0.1) 


9-0(0.7) 


-2(1) 


t < 50 

Myr 
(N=196) 


9.2(0.3) 

11.5(0.3) 

8.6(0.2) 


-26.7(0.4) 


3.24(0.22) 


-0.19(0.07) 


-1-4(0.2) 


9-5(0.7) 


0(1) 


t > 50 
Myr 
(N=178) 


9.4(0.3) 

10.0(0.2) 

7.8(0.2) 


-25.5(0.5) 


5.57(0.37) 


-0.89(0.15) 


-3.5(0.2) 


5-6(0.3) 


-1(1) 
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Table 2: Kinematic parameters of young (< 50 Myr) OSCs in the arms 



Arm 


^0,^0,^0, 
km/s 


km/s/kpc 


^0, 

km/s/kpc^ 


km/s/kpc^ 


K, 
km/s/kpc 


kpc 


^0) 

deg. 


Car-Sag 
and Perseus 
(N=92) 


9.0(0.7) 
9.6(1.4) 
9.5(0.4) 


-27.3(0.6) 


4.04(0.41) 


-0.37(0.16) 


-1.8(0.4) 


8.3(1.0) 


1(1) 


Orion 
(N=93) 


8.7(1.0) 
14.2(0.4) 
9.5(0.3) 


-31.6(1.0) 


4.46(0.84) 


0.23(0.38) 


-0.5(0.6) 


5.4(0.9) 


-3(2) 


Car-Sag 
and Perseus 

(N=92) 


6.5(0.5) 
8.7(1.1) 

9.1(0.4) 


-25.8(0.5) 


4.52(0.37) 


-0.70(0.17) 


-2.6(0.4) 


7.1(0.6) 


1(1) 


Orion 
(N=93) 


8.7(0.5) 
14.0(0.4) 
8.3(0.3) 


-31.7(0.6) 


5.14(0.65) 


-0.20(0.34) 


-0.5(0.4) 


4.9(0.5) 


-3(1) 


Car-Sag 

and Perseus 
(N=92) 


6.5(0.6) 
8.7(1.0) 


-25.9(0.5) 


4.55(0.38) 


-0.69(0.17) 


-2.6(0.4) 


7.1(0.7) 


1(1) 
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Fig. 1. Distribution of OSCs (a) younger than 50 Myr and (b) older than 50 Myr in the 
Galactic XY plane. 
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r,pc 

Fig. 2. Mean errors ay^ (1), cfvi (3), and cxy^ (2) versus heliocentric distance. 




Fig. 3. Galactic rotation curve from clusters within 2.5 kpc (1) and 5 kpc (2) of the Sun. 
The thin lines mark the boundaries of the la confidence intervals; the vertical line indicates 
the position of i?o = 7.4 kpc. 
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Fig. 4. Linear contraction velocity versus (a) heliocentric distance and (b) mean sample age. 



14 



